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ABSTRACT 

Recently it has been reported that there may be a discrepancy between 
big bang nucleosynthesis theory and observations (BBN crisis). We show that 
BBN predictions agree with the primordial abundances of light elements, ^He, 
D, ^He and ^Li inferred from the observational data if the electron neutrino 
has a net chemical potential S,ue due to lepton asymmetry. We estimate that 

= 0.043l[l:[|^[i(95% C.L.) and n.h'' = 0.015 (95% C.L.). 



Subject headings: abundances — early universe — elementary particles — 
nucleosynthesis 



- 2 - 



1. Introduction 

Standard big bang cosmology gives us a simple and natural picture of our universe 
since it explains the origin of the cosmic microwave background and the abundances of the 
light elements, '^He, D, ^He and ^Li. The standard big bang nucleosynthesis (SBBN) has 
a single free parameter, i.e. the baryon to photon number ratio rj = riB/n^ where n-s and 
n.y are the number densities of baryons and photons. Comparing SBBN predictions with 
the abundances of light elements inferred from observational data, we obtain the baryon 
density of the universe, fib ~ 0.01 ([Walker et al. 199l| , [Kernan fc Krauss 1994|) . 



However, recently Hata et al. (1995) pointed out that there is a discrepancy between 
theory and observation by using detailed statistical analysis for both observational data 
and theoretical prediction. This discrepancy is called BBN crisis. Of course if we adopt 
the large uncertainties of the D data observed in Ly-a clouds ( ^ongaila et al. 1994 , Rugers 



Hogan 1996a| , [L996bD , BBN predictions agree with the observed abundances ( |Ulive & 



rhomas 1996| , [Copi, Schramm fc 'lurner 1995| ) . However if we adopt the chemical evolution 



model of D and '^He or low D data observed in Ly-a clouds ( [lytler. Fan fc Buries 19"96| , 
Buries fc Tytler 1996|) , the crisis still exists. In this paper we discuss the problem by 



taking the discrepancy between the theory and observations seriously. One way to solve the 
discrepancy might be to adopt some modifications of standard physics used in SBBN, for 
example, a large chemical potential for a neutrino or an O(lOMeV) mass ( Kawasaki et al. 



|1994|) . It is also pointed out that an exotic massive decaying particle may destroy D and 
^He after the BBN epoch ( [Holtman, Kawasaki fc Moroi 1996| ). 

Since there is baryon asymmetry in the universe, it is quite natural to consider the 
possibility that there is also lepton asymmetry.0 If there exists a lepton asymmetry for 



^In the standard model during electroweak phase transition the large lepton asymmetry is converted into 



baryon asymmetry ( Shaposhnikov 1992). Then the lepton number density is roughly the same as baryon 
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neutrinos, they will have finite chemical potentials. One can reduce the abundance of ^He 
if an electron neutrino has a large chemical potential. This neutrino degeneracy was first 
studied by ([Wagoner et al. 1967|) . So far many people have imposed the constraints on 
the neutrino chemical potentials for z/g ( [Terasawa fc Sato 1988] , |1984| , |1985| , [Reeves 1972 



Beaudet fc Goret 1976[ , [Sherrer 1983[ [Kim fc Lee 1995[) , for z/g and i/^ ([Yahil fc Beaudet 



1976[ , Peaudet fc Yahil 1977[ ) and for three neutrinos Ve-, ([Kang fc Steigman 1992 



Olive et al. 199 1[) . However a precise statistical analysis was not made in the previous work 
and the observational data have been updated since that time. Therefore in this paper, 
we study BBN with the effects of the neutrino degeneracy using Monte Carlo simulation 
( [Smith, Kawano fc IVlalaney 1993 ) and make a likelihood analysis ( Hata et al. 1995 ) using 
the most recent data. 



2. BBN with Neutrino Degeneracy and Likelihood Analysis 

Theoretical predictions should be compared with the constraints of observational data. 
Hereafter we use the notation Yp ( ''He mass fraction), = (D/H)p, y^p = (^He/H)p, 
= ('^Li/H)p , (number fraction relative to H respectively) where "p" denotes the primordial. 
For ^He we adopt the recent analysis of observations in low metallicity extragalactic HII 
regions, from which primordial abundance is given by ( [Olive, Skillman fc Steigman 199^ ) 

Yp = 0.234 ± 0.002(stat) ± 0.005(syst) (1) 

It is believed that the abundances of ^Li observed in the population H metal poor 
halo stars in our galaxy represent the primordial values. We adopt the most recent data 

number density ub ^ 10^ ^"^n^. Thus the total lepton number asymmetry is smaU. However lepton number 
asymmetry for each species of neutrino may be large. 



( IBonifacio fc iVlolaro 199BD P1: 



yrp = (1.73 ± 0.05(stat) ± 0.20(syst)) x 10 
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(2) 



Finally we discuss the constraints of D and ^He. The abundance of D observed in the solar 
neighborhood and inter-stellar matter give us lower limits to the primordial abundance, 
which leads to the constraints: 



?/20 < y2pRx, 



y2,ism < ^20, 



(3) 
(4) 



where Rx = Xp/X„o^ and X is the mass fraction of the hydrogen ( the subscripts "ism", 
"©"and "raoif" denote inter stellar matter, pre-solar, the present universe, respectively ). 
We also adopt a constraint on D and ^He by using the chemical evolution model which 
gives a relation between the primordial and pre-solar abundances of D and ^He ( |Steigman 
k Tosi 19951): 



-.2.+ --lh/3. 



9?. 



{yip + y3py2p)Rx < o. 



(5) 



where is the survival factor of ^He, and is estimated to be 0.25 - 0.50 ( Pearborn" 



Schramm fc Steigman 1986|) P|. The observational abundances with la error are given by 
( peiss 1993| , [Linsky et al. 1993t ^teigman fc Tosi 19951) 



1/20 
1/30 

Rx 

y2,ism 



(2.57 ±0.92) X 10"^ 
(1.52 ±0.34) X 10~^ 
1.1 ±0.05, 
(1.6 ±0.2) X 10~^ 



(6) 
(7) 
(8) 
(9) 



^Here systematic errors due to cosmic spallation and depletion in stars are not included. However, even 
if we include such errors the conclusion in this paper is not changed. 

•^The survival factor may be greater than previously thought due to the uncertainty in the estimation 



of the galactic chemical evolution ( Schramm 1996 ), however we take the conservative values in this paper. 



- 5 - 



Here we make a comment on the recent measurements in Ly-a absorption lines of 
QSOs. It was expected that they would provide D abundances close to the primordial 
value. However two groups have reported inconsistent results on the measurement of D/H. 
{y2,Ly-a ~ (1-9 - 2.5) X 10"^ QSongaila et al. 1994| , [Rugers fc Hogan 1996a| , [T996E|) . and 
y2,Ly-a ~ 2.5 X 10"^ ( rTytlcr, Fan fc Buries 199^ , [Buries fc Tytler 199q ). The lower value 
is the same order as that obtained by the measurements in the solar system. The higher 
value leads to better agreement between theory and observation ([Olive fc Thomas 1996| ). 
At present it is premature to judge which result is correct. Therefore we adopt neither of 
these values in this paper. 

The lepton asymmetry L = J2 Li is defined in analogy with the baryon asymmetry as 



L, = 



ni - nj 



(10) 



where is photon number, rii is number density of lepton I {I = e, fi, r, z/g, z/^, Ur) and I 
denotes the anti-particle of /. If the lepton / is in thermal equilibrium with temperature Ti 
much higher than its mass, rii and nj are given by 



ni 



91] 



X 



9i: 



It 

27r^ J 1 + exp(x — ^/ 



-dx, 



(11) 
(12) 



'2ti^j l + exp(x + ^j; 
where = fii/Ti denotes the degeneracy parameter of each lepton, fii is the chemical 
potential and gi is 1 for neutrinos and 2 for charged leptons. Then the energy density of 
each type of lepton is given by 



Pi = 9i 



X 



+ 



X 



dx. 



'27r^ J yi + exp(x — ^/) l + exp(x + ^/;y 
Since we only need |^;| ^ 1 to solve the discrepancy in BBN, the energy densities are 
expanded in as 

7 /27r2\ . 



(13) 



Pi = 9i 



8 V 30 



^ 30 f^iY 15 f^iY 



(14) 
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Since pi is even function of and monotonically increases, a non-zero speeds up the 
cosmic expansion compared with the non-degenerate case, which leads to the earher n/p 
freeze out and more ^He production. 

The lepton to photon ratio is also expanded in as 



The pre- factor of the above equation is 0(1). Thus the lepton to photon ratio is roughly 
equal to the sum of the degeneracy parameters At the epoch of the nucleosynthesis 
(T ^ 1 MeV), muons and tau leptons have already disappeared through annihilation and 
decay processes. Thus the degeneracy parameters and are effectively zero. Since 
the neutrality of the universe tells us {rie — n^jn^ ~ nsln^ = V ~ 10^^° , is also 
negligible. Therefore we need pay attention to only the lepton number asymmetry of the 
neutrinos. Among degeneracy parameters of three species of neutrinos { i/g , i/f^ , Vt ) , ^i,^ is 
the most important since the electron neutrino degeneracy modifies the neutron to proton 
ratio at freeze-out. Using the condition for chemical equilibrium, the neutron-proton ratio 
{n/p)f at freeze-out is given by {j) j, — ^w{~Q /Tf — where Tf denotes the freeze-out 
temperature and Q — 1.29MeV. Thus if is positive, neutron fraction is suppressed and 
hence the production of *^He is reduced. This effect is more significant than the speed-up 
effect due to increase of p^,^. Thus we expect that the abundance of ^He decreases if > 0. 
On the other hand the degeneracy parameters and have only the speed-up effect 
which increases the n/p ratio. Therefore large mu- or tau- neutrinos degeneracies or 

^ ^i/e) might compensate the decrease of the n/p due to electron neutrino degeneracy. 
However it seems unnatural for ^i,^ or ^i,^ to be much larger than ^j,^ , and it is easy to 
reinterpret the results for — — in the case of large ^j,^ or In this paper we only 
study the big bang nucleosynthesis with the electron neutrino degeneracy. 

We use Monte Carlo simulation and the hkelihood method to obtain the best fit values 
for ^i,^ and rj and estimate their errors. We first compute the light element abundances 




(15) 
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and their uncertainties, using the Monte Carlo method. The uncertainties come from 
statistical errors in measurements of the neutron lifetime ( [Particle Data Group 199BD and 



nuclear reaction rates ([Smith, Kawano fc Malaney 1993]) . We assume that the distribution 
of reaction rates is gaussian. Then the distributions of the calculated abundances of light 
elements also become approximately gaussian. The theoretical likelihood function of the 
light element i is given by 

/ 



LT{yi,U.r]) = -exp 



(16) 



where the mean value /ij and the variance af depend on the model parameters, r] and 
and i = 2, 3, 4, 7 denote the D, ^He, ''He, and ^Li, respectively. 

Let us estimate the observational likelihood functions. For observational data of "^He 
and ''Li we have the statistical errors agtat and the systematic errors (Jgyst ■ Here we treat 
the statistical uncertainty as a gaussian distribution, and the systematic one as a top hat or 
flat distribution. Namely for y"^'^ = ^f"^ ± a°Jl^^ ± ob^^ , the likelihood functions are given by 



sys- 



2 (afyl+ + a 



obs 
sys- 



,(17) 



where erf(x) is the error function. The full likelihood function of "^He and ^Li are obtained 
by convolving the theoretical and observational likelihood functions: 

LiiU, V)=l dy.Lf{y,, v) x Lf{y,). (18) 
Furthermore we regard the likelihood function as a probability distribution function. 

For D and ^He , we take a more complicated procedure. The statistical error is treated 
as a gaussian distribution. Since we use the observed pre-solar abundances of D and ^He, 
we should perform the integration over regions V in which the solar and galactic evolution 
constraints, Eqs.(|D, (^ and (^ are satisfied for a fixed ^He survival fraction. We take 
= 0.25 here ( [Hata et al. 1995| , [Holtman, Kawasaki fc Moroi 19961) . Thus the likelihood 
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function for D and ^He is given by 

^23{^u,,v) = / dy2ismdy2Qdy3QdRxdyf dyf'L2ism{y2ism) x ^20(1/20) x L^Q^y^^) 
Jv 

xLnARx) X Lf{yl\f,f{U,v),4\U,v)) x Lf{yf,fif{U,v),4\U,v))- (19) 
Then the total hkehhood function is given by 

LtotiU,v) = L^iU^v) X L23iU,v) X Lr{^^^,r]). (20) 
The best fit region is obtained by the method. 

The three hkehhood functions Li,L23 and Lj are shown in Fig.l for SBBN, i.e. C,ue 
= 0. It is seen that these functions do not agree with each other. In particular there is a 
significant difference between L4 and L23, which leads to the BBN crisis pointed out by 



Hata et al. 1995| . The positive electron neutrino degeneracy reduces abundance of the ^He 
as explained before. In fact in Fig.l (A) it is seen that the Yp likelihood function moves to 
the right in rj for C,ue = 0.05. On the other hand the likelihood functions L23 and L7 are 
not changed very much since they are almost independent of the neutron abundance at 
freeze-out. Then the four likelihood functions agree with each other and the discrepancy in 
SBBN is solved (Fig. 2). In Fig.3 we show the contours of the confidence levels in the rj-^^^ 
plane. The best fit values for rj and are estimated as 

r^io = 4.0 IJI (95% C.L.) , (21) 

U = (4.3 ^to) X 10"2 (95% C.L.) , (22) 

where rjio = r] x 10^°. Notice that standard BBN ( C,ue = 0) is outside of the 95% C.L. 
contour. Though we adopt a lower here, in the case of a higher survival factor {g^ = 0.5) 
the result is not changed very much, see Fig.2. 
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3. Conclusion 

In this paper we have investigated BBN with lepton number asymmetry. Assuming 
that the electron neutrino has a significant chemical potential due to lepton asymmetry, it 
has been shown that the observational data agree with the theoretical predictions and the 
BBN crisis is solved. We have estimated the lepton number density of electron neutrinos and 
the baryon number density. The baryon density parameter is given by fi^/i^ = 0.015 -o'xm 
(95 % C,L.), where h is Hq = 100 h km/s/Mpc. The estimated chemical potential of z/g is 
about 10~^ eV which is much smaller than experiments can detect (~ 1 eV, [Konstantinov fc 
[Kocharov 1964| ). In other words, BBN gives the most stringent constraint on the chemical 
potential of Ug. 

We thank P. Kernan for helpful advice about BBN code, and T.Totani for valuable 
discussions about statistical method. This work has been supported in part by the 
Grant-in-Aid for COE Research (07CE2002) of the Ministry of Education, Science, and 
Culture in Japan , and in part by Grant-in-Aid for Science Research Fund of the Ministry 
of Education, Science, and Culture No. 08740194. 
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Fig. 1. — (A) Likelihood functions of Yp ("^He mass fraction), (B) likelihood functions of 
D & ^He for — 0.25, 0.50, (C) likelihood functions of '^Li. The solid line represents for 
^i,^ — 0.00 (SBBN) and the dashed line represents for — 0.05. For ^He with an electron 
neutrino degeneracy, the peak of the function moves to the right in 77 axis. The intensity is 
an arbitrary scale. 

Fig. 2. — The contours of the confidence levels in ( abundance — rj ) plane for each light 
elements ^He, D and ^Li for ^^.^ — 0.05. In the contours the solid line denotes 68%C.L. and 
the dotted fine denotes 95%C.L.. Theoretical predictions are also shown with theoretical 
uncertainties (1 a) obtained by Monte Carlo method (for gs = 0.25). 

Fig. 3. — The contours of the confidence levels in the parameter space, r) - for g^ — 0.25, 
0.50. The dashed (sohd) contour corresponds to 68% C.L. (95% C.L.). Standard BBN ( 
= 0.00) is outside of the 95% C.L.. 



